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Angle-resolved ultraviolet photoemission spectra are in- 
terpreted by combining the energetics and spatial properties 
of the contributing states. One-step calculations are in ex- 
cellent agreement with new azimuthal experimental data for 
GaAs(llO). Strong variations caused by the dispersion of the 
surface bands permit an accurate mapping of the electronic 
structure. The delocalization of the valence states is discussed 
analogous to photoelectron diffraction. The spatial origin of 
the electrons is determined, and found to be strongly energy 
dependent, with uv excitation probing the bonding region. 

PACS numbers: 79.60.-i, 61.14.Dc, 73. 20. At 



Angle-resolved ultraviolet photoemission spectroscopy 
(ARUPS) is probably the most powerful single exper- 
imental technique for studying the valence electronic 
structure of solids. However, until recently, ARUPS has 
been primarily limited to investigating the positions in 
energy of valence bands along a few high-symmetry di- 
rections in reciprocal space. That is, the intensity of the 
photoemission peaks was usually not analyzed quantita- 
tively, and most of the information in the full hemisphere 
above the surface was lost. One reason for this limita- 
tion in ARUPS studies is the lack of any simple rules for 
explaining such valence spectra beyond those that have 
been found useful for mapping bands in energy. The 
most often applied model in ARUPS is that of direct 
(wave-vector-conserving) transitions between bulk bands 
[Ql, with the final state often being simplified to a plane 
wave Beyond this, free-electron final states with 

atomic-like optical transitions jsj and final-state scatter- 
ing of electrons emerging from a localized core orbital ^] 
have been used to better understand the resulting angu- 
lar distributions in photoemission. The interest in such 
angular distributions of intensity has recently been stim- 
ulated by the measurement of full-hemisphere intensity 
maps in ARUPS [p|-p^, with the data being analyzed so 
far using plane- wave final states , single scattering , 
bulk bands H , and atomic-like transitions |^ , as well as 

The most accurate descrip- 



the one-step model |11 



elements, full multiple scattering, the explicit presence 
of the surface potential, and the resulting more complex 
initial and final states. However, there has to date been 
no systematic treatment of the angular distribution of 
intensity within such a model. 

In this Letter, we quantitatively investigate the influ- 
ence of the initial state and its charge distribution on 
the angular distribution of intensity in ARUPS, using 
the one-step model. We demonstrate that the combined 
consideration of both energy positions and intensity pat- 
terns will be necessary for the most useful interpretation 
of angle-scanning data, and point out the additional in- 
formation that can be derived in this way. For our calcu- 
lations and measurements, azimuthal scans are chosen, 
which allow a high accuracy in the visualization of the 
angular distribution. A generalization of the x-ray pho- 
toelectron diffraction (XPD) picture is found to lead to 
a more unified view of the electronic structure in both 
direct and reciprocal space. The GaAs(llO) surface is 
chosen as a well-understood test case |l^^-[lq| for which 
it will become evident that the photoemission patterns 
reveal not only the surface density of states (SDOS) but 
also provide insight into the charge density of the bonds. 
This contributes to a common understanding of ARUPS 
and XPD angular distributions. 

Our interpretation of photoelectron emission proceeds 
via the construction of the photoelectron state ^ with 
energy which can be written at the detector lo- 

cated in the direction of the polar and azimuthal (ip) 
emission angles as 



DGiEf,r.,i9,ip) (A-p) I*,, 



(1) 



tion of valence ARUPS involves calculations within this 
one-step model and includes the precise optical matrix 



where D is the SDOS, G is the propagator, A is the vec- 
tor potential, and p is the momentum operator Q . The 
initial state, depending on the energy Ein and k|| , is split 
into the factor of the SDOS D and the wave function "ifin- 
This formula is used for the following discussion of the 
electronic and spatial structure. Other influences like se- 
lection rules, the density of the final state, or resonances 
of direct volume transitions are not explicitly apparent, 
but they are all included in Eq. (0) and in a one-step 
calculation for which Eq. (|l|) is transformed to a golden 
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rule formula. There is principally no difference in the 
validity of applying the one-step model to both ARUPS 
or XPD. In Eq. (|^), D{Ein,'k\\) contributes to the band 
mapping, and, as noted above, is often the only thing 
considered in most analyses to derive experimental band 
structures. G describes the electron scattering from the 
point of excitation through the crystal to the detector. 
In XPD this scattered electron can be simply thought 
as being emitted from a localized initial state; here the 
general case arises to treat emission from delocalized va- 
lence states. The source of the scattered electron is no 
longer the vicinity of a core but the entire initial state 
wave function. The wave amplitude (r|^'i„) in Eq. (|^) 
gives the local emissivity, depending on the spatial po- 
sition, initial energy, parallel momentum, and, together 
with the radiation polarization, the angular momentum 
composition. 

As a test of this interpretation, we investigate 
azimuthal distributions of the photocurrent above 
GaAs(llO) in terms of the electronic and geometric struc- 
ture of the participating states. The calculation of the 
photocurrent within the one-step model proceeds along 
the lines described elsewhere The spectra were 

measured with unpolarized HbIq, radiation {hv — 21.22 
eV) with the energy resolutions set to 130 meV and, for 
testing purposes, to 35 meV, which gave negligible differ- 
ences. The scans were taken with an angle resolution of 
0.25° by use of a 180° spherical electron analyzer, which 
is movable around two independent axes. These degrees 
of freedom permit taking angle scans without moving the 
sample and allow studying the influence of the incoming 
light and separating such effects from other processes. 

Two cases are studied here: At —0.75 eV there is a 
dangling-bond state (A5), that is a surface state lying 
outside the projected bulk bands, and at —4.0 eV there 
are two resonances, the Ga back-bond state (C2) and the 
A3 state. To get a broad view off the high symmetry 
points, we consider a circle in the Brillouin zone with a 
radius of fc|| = 0.6 A~^, corresponding to the measure- 
ment angles at both energies. The surface band structure 
along this circle is shown in Fig. ^together with the pro- 
jected bulk band structure. The energy contour plot in 
the inset of Fig. |l| exhibits the symmetry of the A5 state; 
this surface has only a mirror plane, but as an additional 
symmetry, an inversion symmetry occurs in the electronic 
structure due to Kramer's degeneracy (-E(k) = k)). 

For a detailed discussion of the electronic and geomet- 
ric influences on the photocurrent, the effect of the inci- 
dent light on these spectra must be known first. With a 
fixed polar angle of incidence and the use of unpolarized 
light, an asymmetry remains only due to the azimuthal 
direction iphu of the light. Here, enhancements and az- 
imuthal peak shifts are weak, as shown in Figs. || (a)-(c). 
These influences can be discriminated from the strong 
differences in the current between positive and negative 
y directions. 



For this A5 state, a comparison of the photocurrents 
in Figs. I (a)-(d) with the SDOS in Fig. | (e) shows that 
the SDOS controls the overall structure in the current. 
The number and the azimuthal positions of the lobes are 
the same for the SDOS and the current, but differences 
occur in intensities and in smaller structures. There are 
4 broad lobes in the SDOS and the current, but instead 
of 8 maxima as for the SDOS on a fine scale, the cur- 
rent displays only 4 maxima. Figs. || (a) and || (b) for 
the C2-I-A3 state show the same overall correspondence 
and differences in intensities. The SDOS used here is the 
partial density of states of the orbitals from the upper- 
most atomic layer. The power of the SDOS for accurate 
investigations of the energetic structure is demonstrated 
with the emission from the A5 state. Here, this excellent 
agreement between experimental and theoretical currents 
could only be achieved with a small shift in the corre- 
sponding binding energies by 0.15 eV, but the effect of 
this energy shift is remarkable, as shown in Fig. ^ (d). 
The main reason for this difference should be attributed 
to the calculated surface band structure, though the in- 
accuracy is still in the usual range of common theoretical 
uncertainty [ p^ . The strong variation of the shape in 
Fig. I (d) reflects the changes in the SDOS in Fig. | (e). 
This is obviously caused by the flat dispersion of the band 
which is shown in Fig. |l|. Thus the number, position and 
intensity of the lobes in the current and SDOS depend 
very sensitively on details of the band structure. Small 
shifts in the energy cause not only shifts in the accom- 
panying peaks, but can give rise to strong changes of the 
entire pattern. This allows an accurate determination of 
the band structure by comparing experimental with the- 
oretical scans. The usual simpler approach of identifying 
the maxima in the current as band positions, would fail 
here to give even the correct number: For the A5 state 
the band is hit 8 times (cf. inset in Fig. 0), which corre- 
sponds to the 8 maxima in the SDOS, but the current in 
Fig. ^ has only 4 maxima, and one would overlook half 
of the band positions. 

In contrast to the electronic structure, inversion sym- 
metry is lacking in the photocurrent. To understand this, 
we have to consider in addition to the SDOS also the wave 
function (r|5'i„) of the initial state, which are calculated 
in a LCAO basis. For the dangling-bond A5 band, the 
asymmetry between the intensities in positive and nega- 
tive y directions corresponds to the direction of the bond. 
The dangling bond points along the negative y direction 
and into the vacuum in the z direction, as shown in Fig. |^ 
(f) by the contour plot of a As-charge density. Since the 
density is located above the uppermost As atoms, there is 
only weak potential scattering for a major part of the ex- 
cited electrons. Therefore, the SDOS and the initial state 
wave function dominate the angular distribution, and a 
free final plane wave may be a sufficient description. The 
dotted current in Fig. |^ (d) shows how well this approx- 
imation does in reproducing the general asymmetry. Be- 
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cause a Fourier transformation of a spherical harmonic 
reproduces the same spherical harmonic and because of 
the simple spatial and angular momentum structure of 
the dangling-bond state, the asymmetry of the current 
can be identified here with that of the charge distribu- 
tion. Such a connection has already been assumed in 
a former ARUPS analysis and it is quantitatively 
proven here. As a result, the photoemission intensities 
reflect directly the energetic and spatial density of the 
initial A5 state. However, this simplicity is not a general 
rule, as we will illustrate for our C2-I-A3 case at —4.0 eV. 

To achieve further insight into the influence of the ini- 
tial state, the importance of different parts of the emit- 
ting volume is investigated in a novel way. Apart from 
the charge density maxima in the bonding region the ini- 
tial wave functions also have charge peaks closer to the 
cores, as shown in Fig. ^ (f) and Fig. ^ (f). In the present 
case, three areas are separated by almost spherical nodal 
surfaces. The separate contributions to the current from 
each of these regions are calculated by setting the initial 
wave function to zero in the other two areas. Contribu- 
tions for the emissions from —4.0 eV are shown in Figs. ^ 
(d)-l-(e). From the localized middle region arises less than 
1% of the total current for both studied binding energies. 
Contributions from the innermost area are even smaller 
and completely negligible. With a free flnal state for 
simplicity we have repeated these calculations for kinetic 
energies up to 1500 eV. Below 100 eV only contributions 
from the outer volume are notable, but at 200 eV the 
current from the middle area becomes appreciable. This 
middle area dominates above 500 eV, and constitutes 
the current above 1000 eV. The spatial scale detected 
by photoelectrons may be discussed in close analogy to 
the energy dependence of common scattering processes 
where high wave vector Fourier components probe the 
more rapidly varying wave functions near the core. Usu- 
ally with increasing energies shrinking vicinities of the 
cores are studied. 

At —4.0 eV, neither the number nor the position of 
the lobes nor the asymmetry in the intensities is repro- 
duced with a final plane wave, as shown in Fig. ^ (c). 
Only the inclusion of scattering provides the asymmetry 
in the intensities, as the excellent agreement in Fig. || (a) 
between experiment and theory shows. Therefore, the 
simple argument with the Fourier transformation as in 
the A5 case does not work. We trace the origin of the 
intensities into the details of the calculation. Although 
the A3 state extends over 4 atomic layers, 80% of the 
electrons are excited from the density around the up- 
permost As atom and from the bonds towards the next 
Ga atoms. The inversion asymmetry of the main peak 
at 15° is caused by these bonds, whose unsymmetrical 
charge distribution is shown in Fig. |^ (f). This lack of 
the band-structure symmetry in the intensities generally 
seems to give a strong hint as to the asymmetry of the 
charge density. 



The dependence of the spatial distribution on kinetic 
energy gives an additional feature of ARUPS from va- 
lence states and its relation to XPD. The importance of 
smaller spatial scales with increasing energy might be un- 
derstood as a localization of an effective emitting source. 
At high energies this coincides with the observation that 
the XPD pattern from valence bands are nearly identical 
to those from the localized core states ||l^]. The com- 
mon spatial origin of the excited electrons leads to sim- 
ilar currents since in the one-step model the flnal states 
are exactly the same. It should be noted that at these 
high energies the flnal state scattering is known to dom- 
inate over further details of the source wave 1 18 1^ , and 



especially initial state interferences vanish on the aver- 
age by the flnite resolution Contrarily, at low ki- 
netic energies valence states may contribute from regions 
where the localized core states vanish. These delocalized 
regions are the most important for both studied cases. 
Whereas with the effective localization at the core the 
XPD patterns reflect the geometry, ARUPS probes the 
bonding region with an intensity distribution showing in- 
formation about the bonds. This is a new aspect for the 
interpretation of ARUPS data. 

We have presented a joint treatment of the spectral 
and spatial features of ARUPS from the model surface 
GaAs(llO). The azimuthal scans reflect in number and 
angular positions of the lobes the electronic structure 
given by the SDOS. By comparing measured and calcu- 
lated currents, huge changes induced by dispersion allow 
an accurate determination of the SDOS. Even for unpo- 
larized light there are strong additional intensity modu- 
lations, which are connected to the initial wave function. 
For an interpretation we adopted the XPD picture gen- 
ralized to delocalized valence states. The amplitude of 
the initial state appears as the local emissivity for the 
spatially distributed source of the electrons to be scat- 
tered. Photoelectron spectroscopy is sensitive to different 
spatial parts of the initial state, depending on the kinetic 
energy. Ultraviolet photoemission detects the wave func- 
tion in the bonding region outside the core, and the origin 
of photoelectrons are traced for the flrst time into single 
bonds. This opens new possibilities in the application 
and interpretation of ARUPS. 
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FIG. 1. Calculated surface band structure with the peaks 

(circles) of the SDOS along the circular path for = 0.6 
A^^. Strong peaks are dark. The projected bulk band struc- 
ture is also shown (shaded). The inset shows the As band in 
a contour plot with equidistant energy levels and the level for 
—0.75 eV (dark). The band disperses in the Brillouin zone 
from r down to a local minimum in X' and its global mini- 
mum in X. The dotted circle describes the ky path studied 
here. All energies are referred to the valence band ma^ximum. 



(c) 98°|'y, k 




FIG. 2. Photocurrent at = 0.6 A"^ (i? = 17.3°), SDOS, 
and charge density of the A5 band: (a) to (c) current mea- 
sured at —0.9 eV (black) and calculated at —0.75 eV (grey) 
for light incident at azimuthal angles of 98° , 228° , and 8° , re- 
spectively; (d) calculated current for 98° incidence at —0.75 
eV (grey), at —0.9 eV (solid), and at —0.75 eV with a free 
final state (dotted); (e) SDOS at -0.75 eV (grey) and at -0.9 
eV (black); (f) contour plot of the charge density in a yz plane 
at —0.75 eV for ky belonging to the lower left lobe as indi- 
cated by the dashed lines in (d) and (e). In the polar plots 
(b) to (e) , the radius represents the intensity of the current 
and the value of the SDOS. Accordingly, the azimuthal angle 
is referred to the emission direction and to the ky vector. The 
orientation of the axes is indicated in (c). 
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FIG. 3. Photocurrent at -4.0 eV and = 0.6 A'^ 
(i? = 19.4°), SDOS and charge density for C2+A3 states: 
(a) measured (black) and calculated (grey) current for light 
incident from 88°; (b) SDOS; (c) current with a free-electron 
final state; (d) + (e) partial currents emerging from the mid- 
dle (d), and the outer (e) area for light incident from 90°; (f) 
contour plot of the charge density in a plane through the the 
uppermost As and Ga atoms, at —4.0 cV, and for ky belong- 
ing to the upper right lobe as indicated by the dashed lines 
in (a) and (b). The current in (d) is magnified by 100 rela- 
tive to that in (e). The middle region is the spherical volume 
between radii of 0.08 A and 0.29 A around the As cores and 
0.09 A and 0.3 A around Ga as depicted in (f). The outer 
region is the space outside these spheres. 
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